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ABSTRACT 

The microquasar GRO J1655-40 has a black hole with spin angular momentum ap- 
parently misaligned to the orbital plane of its companion star. We analytically model 
the system with a steady state disc warped by Lense-Thirring precession and find the 
timescale for the alignment of the black hole with the binary orbit. We make detailed 
stellar evolution models so as to estimate the accretion rate and the lifetime of the 
system in this state. The secondary can be evolving at the end of the main sequence 
or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in 
the latter case but we find that, in both cases, the lifetime of the mass transfer state is 
at most a few times the alignment timescale. The fact that the black hole has not yet 
aligned with the orbital plane is therefore consistent with either model. We conclude 
that the system may or may not have been counter-aligned after its supernova kick 
but that it is most likely to be close to alignment rather than counteralignment now. 
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1 INTRODUCTION 

Microqu asars are X-ray binaries wh ich have relativistic ra- 
dio jets (|Mirabel fc R odriguez 1999). They consist of a com- 
pact object such as a black hole or a neutron star which 
is accreting matter from a companion star. GRO J1655-40 
is one of only fifteen observed microquasars in our galaxy 
(Parcdcs 2005) and has a black hole of mass Mi = 6.3 ± 
0.5 Ma and a companion star of mass M2 = 2.4 ± 0.4 Mq 
i|Greene. Bailvn fc Oros3l200lD . 

GRO J 1655-40 shows transient beh aviour in a complex 
form. It was discovered in July 1994 (jZhang et al.l ll994T ) 
through observations of hard X-ray outbursts and then ob- 
served again several times in the following year. Before 1994 
it appeared to be quiescent for some time. The original out- 
bursts were associated with radio flari ng and superluminal 
motio n of the radio plasmoids which led lHiellming fc Rupenl 
[| 19951 ) to estimate the distance to the system to be 3.2 kpc. 

Ilsraelian et al.l (1 19991 ) measured the abundances of oxy- 
gen, magnesium, silicon and sulphur in the atmosphere of 
the companion star to be six times greater than in the Sun. 
The secondary star does not have enough mass to reach 
the internal temperatures required to create these elements. 
A stellar wind or mass transfer from the primary during 
its presupernova evolution could only have provided a small 
fraction of these observed overabundances because the CNO 
cycle would have led to a much larger nitrogen to oxygen ra- 
tio than is observed. So these abundances are interpreted as 
evidence that supernova ejecta have been captured by the 
secondary star, perhaps during fall-back on to the black hole. 
The relative abundances of the contaminating elements led 



Ilsraelian et ail (|l99Sl > to suggest that the black hole's pro- 
genitor was a star with mass 25 — 40 Mq . It is therefore 
likely that the system lost more than half its mass in the 
supernova even if the progenitor of the black hole had had a 
strong stellar wind. In order for the system to have remained 
bound the supernova explosion must have been associated 
with a substantial kick. Such a kick could have altered the 
spin of the black hole leaving it misaligned with the binary 
orbit. Even a small velocity kick can lead to a large off- 
set between the spin of the bla ck hole and the orbital axes 
|Brandt fc Podsiadlowsklll995l ). There is no reason why a 
kick substantial enough to keep the system bound would 
leave the spin aligned with the orbit. 

GRO J 1655-40 has relativistic jets of material leaving 
the system which we assume are generated in the inner parts 
of a black hole accretion disc and so are perpendicular to the 
plane of the disc close to the hole. The combined action of 
the Lense-Thirring effect and the internal viscosity of the 
accretion disc causes the angular momenta of the black hole 
and the inner accretion disc to align. This is known as the 
iBardeen fc Pettersonl (|l975l ) effect. It affects only the inner 
regions of the disc because of the short range of the Lense- 
Thirring effect. The outer parts of the disc tend to remain 
in their original configuration. 

iHiellming fc Rupenl (Il995l ) measured the jet inclination 
to be 85° ± 2° to the line of sight. This differs significantly 
from the binary orbital plane inclination of 70.2° ± 1.9° 
l|Greene. Bailvn fc O rosz 2001). This implies that there is a 
misalignment between the inclination of the black hole and 
the outer parts of the accretion disc and so the accretion 
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Figure 1. We consider the position of the jet fixed relative to the 
line of sight. The cone represents the surface on which the binary 
orbital angular momentum must lie. The smallest possible incli- 
nation between the black hole and the binary angular momentum 
occurs in position a where r) = 15°. The largest inclination be- 
tween the black hole and the binary angular momentum occurs 
in position c where r\ = 165° . 



disc is warped. The inner parts are aligned with the black 
hole spin and the outer parts are aligned with the binary 
orbital plane because of tidal torques. 

In Figure 1 we show the line of sight and the given 
position of the jet. Relative to this jet, the angular momen- 
tum vector of the binary must lie in one of the cones either 
pointing towards or away from us. The smallest possible 
inclination between the black hole and the binary angular 
momentum, rj, occurs in position a where rj = 15° and the 
system is close to alignment. The largest possible inclination 
between the black hole and the binary angular momentum 
occurs in position c where rj — 165° and the system is close 
to counter-alignment. The inclination of the system is in the 
ran ge 15° < rj < 165° . 

King etaP (|2005T ) find that for a black hole that is mas- 
sive compared to a misaligned disc, then the disc aligns 
with the black hole. If the disc is massive compared with 
the black hole then the black hole aligns with the disc. In 
GRO J1655-40 the disc is much less massive than the black 
hole. However because the disc is in, and is tidally linked 
to the orbitial motion of, a binary system, what matters is 
the amount of angular momentum in the binary orbit, not 
just in the disc. The orbital angular momentum is several 
hundred times that of the black hole so we expect the hole 
to t end to event uall y alig n with the orbit. 
Maccarond 



(2002 



the 



model 



of 



iNataraian &: Pringld 1 19981 ) to calculate the timescale 
on which the black hole aligns with the outer disc and 
hence the binary orbit. This assumes that the surface 
density and viscosities are constant. He finds that the ratio 
of alignment timescale to binary lifetime is about 0.3. He 
assumes the binary lifetime is the time to accrete the whole 
star, not the stellar evoluti on time. U si ng th e outburst 
history of GRO J1655-40 iMaccarond (|2002f ) finds the 
alignment timescale to be 8 x 10 8 yr with a companion of 
main-sequence lifetime of 7 x 10 s yr. He concludes that the 
system should be approximately one alignment timescale in 



age and so the black hole should not be fully aligned with 
the disc because it takes more than one alignment timescale 
for full alignment. 

We examine in more detail the stellar evolution of 
the companion star to estimate how long the system 
can remain in this steady state and whether the system 
should be aligned or not. We make use of the analysis of 
iMartin. Pringle fc Tout! l|2007f ) and use power laws in dis- 
tance from the central black hole for the surface density and 
viscosities to find the alignment timescale of the system. We 
consider the possibility that the disc was initially counter- 
rotating. 



2 MISALIGNMENT 

Because the binary is not visually separated we only know 
the inclination of the orbit relative to the line of sight. We 
do not know its projection on to the sky. Given the observed 
inclination of the binary orbit and the jets to the line of sight 
we can find the probability distribution of the angle between 
the two momenta. 

If the position angle, <f>, of the binary relative to the line 
of sight (see Figure[l} is randomly distributed in < cf> < 2ir 
then 



Since 



we deduce 



1 

2^ 



P{rj)drj = P{<t>)d<t> 



(i) 



(2) 



(3) 



Note that each value of rj maps to two values of 4>. If the 
binary angular momentum is pointing towards us we have 
the standard spherical trigonometric cosine formula 



cos rj = cos i cos j + sin i sin j cos < 



(4) 



where i = 85° is the inclination of the jet to the line of sight 
and j — 70° is the inclination of the binary orbit to the line 
of sight. Differentiating we find 



sin 7; 



drj 



sm i sm j sm < 



(•>) 



If the binary orbital angular momentum is pointing away 
from us then we have the relation 

cos rj = cos i cos(-7r — j) + sin i sin(7r — j) cos <f> (6) 



and we find 



sm rj 



drj 



sm i sm(7r — j) sm < 



(7) 



In Figures [2] and [3] we plot the probability distributions 
for the two cases where the binary angular momentum is 
pointing towards and away from us. We assume either is 
equally likely. We see that if the binary points towards us, 
the most likely misalignment angles are the extremes of rj = 
15° (case a in Figure 1) or 155° (case b) and if the binary 
points away then the most likely angles are rj = 25° (case 
d) or 165° (case c). 
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Figure 2. The probability distribution of the angle between the 
jet and the binary orbit, i), when the angular momentum of the 
binary points towards us. 
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Figure 3. The probability distribution of the angle between the 
jet and the binary orbit, r), when the angular momentum of the 
binary points away from us. 



3 WARPED ACCRETION DISCS 

There are two viscosities in the disc, v\ corresponds to the 
azimuthal shear (the viscosity normally associated with ac- 
cretion discs) and V2 corresponds to the vertical shear in 
the disc which smoothes out the twist. The second viscosity 
acts when the disc is non-planar. We assume that we have a 
steady state disc in which i^iE = const and that the surface 



density is a power law 



(8) 



i|Shakura fc Sunvaevlll973h . where R is the spherical radial 
coordinate, Ro is some fixed radius and Eo is a constant. To 
be in steady state, the first viscosity must obey 



v\ = i/10 



R_V 
Ro) 



(9) 



where /3 and vio are constants. We assume that the second 
viscosity also obeys the power law 



I/O 



V20 



(10) 



and V20 is a constant. 

Following the work of iMartin. Pringle fc Tou j (|2007t ). 
we consider a black hole of mass Mi at R = with spin 
angular momentum J. We consider the disc to be made up 
of annuli of width dR and mass 2nT,RdR at radius R from 
the central object of mass Mi with surface density T,{R,i) 
at time t and with specific angular momentum density L = 
(GMii?) 1//2 EZ = LI. The unit vector describing the direction 
of the angular momentum of a disc a nnulus is I = ( l x ,l v , l z ) 
with |i| = 1. We use equation (2.8) of lPringld (|lQQ2T > setting 
dL/dt — and add a term to describe the Lense-Thirring 
precession (the last one) to give 



R dR 



3R d , T . 3 
-dR (VlL) ~2 Vl 



t 1 r,r dl 

L+ r* RL 8R 



+ 



x L 



i? 3 



The Lense-Thirring precession is given by 

2GJ 



(11) 



(12) 



(Ku mar fc Prinridll985h . where the angular momentum of 
the black hole can be expressed in terms of the dimensionless 
spin parameter a such that 



./ 



acM ^ ( — 



(13) 



The black hole spin is also evolving because of the torques 
exerted by the disc so we have 



dJ_ 

dt 



-2vr 



^3 RdR ' 



(14) 



where the integration is done over the surface of the disc. 
IMartin. Pringle fc Tout] (|2007h solve equation (JTTJ to find 
the disc profile and then use equation (|14p to find the 
timescale for alignment of the black hole with the binary 
orbit. For f3 = this is 



^ align 



(0) 



/ acMj 



(15) 



V27rE \v 2 G 

dScheuer fc Feilerll 199fj ). where E and v 2 must be evaluated 
at the warp radius. The warp radius is found by balancing 
the Lense-Thirring precession term with the viscous term 
associated with the second viscosity in equation (|lip to find 



Rw 



vi{R* 



(16) 
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i|Scheuer fc Feilerl 1 19961 ; iMartin. Pringle fe Toutl I2007I I. In 
the next section we consider the size of this radius for GRO 
J1655-40. 

We can eliminate E using the steady state accretion 

rate, 



M = 3ttEza, 



to find 



talign(O) 



3^i / acM\ 



V2M \ v 2 G 
For (3 the alignment timescale is 



*align(/3) = ialign(0)r, 



where 



(1 + /3) 



(2(1+0)) 



72 



1 (wm) cos (i(iTay) 



(17) 



(18) 



(19) 



(20) 



,2(l+/3) / 

ijMartin. Pringle fc Toutl l2007h and i/i, V2 and E must be 
evaluated at the warp radius, -R wa r P - The alignment of the 
system happens exponentially. The time for the system to 
get very cl ose to alignment is actually several alignment 
timescales (|Martin. Pringle fc Toutl 120071 ). In deriving the 
timescale we assume that we can neglect the non-linear term 



LdH/dRf 



-\dl/8R\ . In the following section we calcu- 



late the alignment timescale of GRO J1655-40. 



4 BLACK HOLE INCLINATION 

We know that the inclination of the black hole to the orbital 
plane is in the range 15° < r\ < 165°. We let 77 be the current 
inclination of the black hole to the outer disc. The outer disc 
is in the binary orbital plane. We let 770 be the inclination 
when the system first started transferring mass. If r/o < n/2 
the system started closer to alignment whereas if 770 > n/2 
the system started closer to counter-alignment than align- 
ment. The system is always alignin g so 77 < 770 even if it is 
count er- aligned. We use the work of lMartin. Pringle fc Tout) 
|2007l ) to find that 



sin 77 = sin rjo e 



(21) 



where taccrctc is the time for which the system has been 
steadily transferring mass, if 770 < n/2 and the system 
started nearer to alignment than counter-alignment. How- 
ever, if the system started closer to counter-alignment then 



sin 77 = sin 770 e 

and 770 > 7r/2. We estimate t aC cr 
in Section [5] below. 



(22) 



from binary star models 



4.1 Co-rotating Disc and Black Hole 

If the system started closer to alignment than counter- 
alignment so that 770 < tt/2 then we find the time that the 
accretion must last for, by inverting equation (|21|l . to be 



= Malign log 



sin 770 
sin 77 



(23) 



However if the disc was initially counter-rotating so that 
770 > 7r/2 and now 77 < 7r/2 then we must find the time that 
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Figure 4. The initial inclination 770 of the black hole relative to 
the binary orbit against the time that the disc has been accreting 
in a steady state. The disc has 77 = 15° now. 



it takes to move from 770 to 90° using equation (|22[) and 
then add on the time it takes to move from 90° to 15° using 
equation (|21|) . Then we find 



-talign log (sin 770 sin 77) . 



(24) 



As an example, in Fig. [4] we see how the initial inclination of 
the disc varies with the length of time it has been steadily 
accreting for a disc which currently is at an inclination of 
15°. 



4.2 Counter-rotating Disc and Black Hole 

If the disc is now closer to counter-alignment than alignment 
so that 77 > 7r/2 then, using equation (|22|l . we find the initial 
inclination of the system as a function of the time for which 
it has been accreting to be 



talign log 



sin 77 



(25) 



, sin 770 , 

In Fig. [5] we see how the initial inclination of the disc varies 
with the length of time it has been steadily accreting for for 
a disc which is now at an inclination of 165° . 

As illustrations we examine further these two extreme 
cases. If the inclination of the black hole relative to the bi- 
nary orbit is random then the probability that it formed at 
15° or less is only 0.017. Similarly the probability that it 
formed at 165° or more is only 0.017. It is therefore likely 
that it has undergone some but not excessive alignment and 
we would expect taccrete ~ talign. We note that a completely 
counter-aligned black hole is in an unstable equilibrium so 
it is possible that we now observe the hole counter-aligned 
at an inclination of 165° . However this would have required 
both the unlikely event of forming very close to counteralign- 
ment and the accretion timescale to be small compared with 
the alignment timescale. It is therefore more likely that the 
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Figure 5. The initial inclination tjq of the black hole relative to 
the binary orbit against the time that the disc has been accreting 
in a steady state. The disc has r\ = 165° now. 



hole is now only 15° from alignment. We cannot tell how mis- 
aligned it started but it may well have been counter-aligned 
in the past. 



5 WARPED DISC MODEL OF GRO J1655-40 

In this section we find the viscosities, the warp radius and 
the alignment times cale of GRO J1655-40. The p eriod of 
the system is 2.62 d (|Greene. Bailvn fc Oros3l200ll ). By Ke- 
pler's law the sep aration is Rb = 1-14 x 10 1 2 cm for a total 
mass of 8.7M0. lAbramowicz fc Kluzniakl (|200lh estimate 
that t he dimensionless spin of the black hole a = 0.2 — 0.67, 
while IWagoner. Silbergleit fc Ortega-Rodriguej l|200lf) sug- 
gest that it could be as high as 0.92. 

We take the outer truncation radius of the disc, i? ou t, 
to be the largest possible stable orbit of a test particle in 
this b inary system. We use parameters derived bv lPaczvnskil 
(1 19771 ) for the three body system. In Fig. [S] we plot the trun- 
cation radius against the mass ratio of the two stars. GRO 
J 1655-40 has mass ratio M2/M1 = 2.4/6.3 = 0.38 which 
corresponds to a truncation radius of i? ou t = 0.37 -Rb = 
4.22 x 10 11 cm which is the outer edge of our disc. This 
outer edge is used later to find timescales in the disc. 
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Figure 6. The truncation radius of the disc as a function of the 
ratio of the masses of the two stars. 



(1999) find the viscosity to be 



Vi = 6.40 x 10 ib a? 



R 



Mi 

Mb 



M 



10" 8 M yr- 



10 n cm 



2 -1 
cm g , 



(26) 



where j — 1,2, cij is the usual a-prescription for viscosity 
v — otc a H and M is the accre tion rate on to the black hole. 
We take a\ = 0.2 and a 2 = 2 (jLodato fc Pringlejl2007l ). This 
is equivalent to equations @ and (|10|1 with Ro — 10 11 cm 
and 



^10 = 1.11 x 10 1 



f-V 

V0.2/ 



Mi 



6.3 M 



M 



lO- 8 M yr- 



and 



V20 = 7.03 x 10 1 



15 (012 



2 -1 
cm g 



Mi 
6.3 M, 



(27) 



M 



10- 8 M o yr- 



2 -1 
cm g . 



(28) 



The viscous timescale associated with v\ is the timescale 
on which matter is moved through the disc. For GRO J 1655- 
40 this is 



5.1 Viscosities 

We use the a-prescription to find the viscosity of the disc 
|Shakura fc Sunvaevl Il973l ). In bright X-ray binaries the 
dominant radiative opacity is bound free. Using the steady 
state mass accretion rate, M = 37r^iE, IWiiers fc Pringlg 



1.7 

M 



f-V 

V0.2/ 



lO-^oyr" 1 



Mi 



6.3 M 



R 

Rout 



5/4 



yr. 



(29) 



This is a relatively short timescale and so the disc reaches 
steady state quickly. 
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5.2 Warp Radius 

We find the warp radius to be 
2lj u 



(30) 



^20(-Rwarp/cm)4 

l|Scheuer fc Feileri 1 19961 : iMartin. Pringle fc Toutl l2007t ) so 

that 



-Rwarp — 4.33 X 



41 / Mi 



6.3M 



M 



lO-^oyr" 1 



(31) 



We note that this is not strongly dependent on M and so 
the outbursts of the system do not affect it greatly. We find 
the viscosities and surface density at the warp radius so that 
we can find the alignment timescale of the black hole with 
the disc. 



5.3 Timescale of Alignment 

In Section[3]we derived a formula for the timescale on which 
the black hole aligns with the disc. This assumes that the 
angular momentum of the disc is much greater than that of 
the black hole. In GRO J1655-40 the outer parts of the disc 
are locked to the binary orbital plane by tidal torques. It is 
therefore the angular momentum of the binary orbit, which 
is much greater than that of the black hole, that matters 
because the disc is constantly fed with matter at its outer 
edge.. 

We can now evaluate the timescale of alignment because 
we know the viscosities and the warp radius at which we 
evaluate them. For a steady state disc we have (3 — | (from 
equations and (|26[) ') and so we find 



£align(3/4) =t a l ign (0) 



4\ 1 



T(2/7) 



7 J V2T(5/7) cos (f) 
1.5233*^(0) (32) 



and 



Mi 



6.3 M 



M 



10- 8 M Q yr- 



yr (33) 



so that 



6 . . _ 

Mi 



6.3 M, 



M 

10- 8 M Q yr- 



yr.(34) 



In Fig. [7] we plot the timescale for alignment against the 
spin of the black hole, a, for a range of accretion rates. In 
Fig. [8] we plot the timescale for alignment of the black hole 
against the accretion rate for a range of spins. In the next 
section we consider the time for which the mass transfer can 
take place to see whether the black hole should be aligned 
or not. 
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Figure 7. The timescale of the alignment of the black hole with 
the disc against a for M = 10~ 7 (solid line), 10~ 8 (dotted line) 
and 10~ 9 M Q yr _1 (dashed line) 
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Figure 8. The timescale of the alignment of the black hole with 
the disc against M for a = 0.2 (solid line), 0.5 (dotted line) 
and 0.92 (dashed line). 



6 BINARY EVOLUTION MODELS 

For our equilibrium warped disc model to be valid, the vis- 
cous timescale in the disc (about 1.7 yr as found in sec- 
tion I5.1|l must be much less than the time for which ma- 
terial has been accreting through it. We note that there has 
been time since the 1994 outburst for the disc to regain 
equilibrium. Outbursts have occurred more frequently in re- 
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cent years but the system generally returns to quiescence 
for sufficiently long periods for the disc to regain equilib- 
rium. Because we still see the jet misaligned, the timescale 
for alignment must be long enough relative to the accretion 
timescale. We therefore investigate the evolutionary state of 
the system to determine both the expected mass-transfer 
rate and the time since mass transfer began. 

Early attempts to make detailed evolution models of 
GRO J1655-40 met with the problem that stars at the sec- 
ondary's position in the Hertzsprung-Russell diagram ought 
to be in the Hertzsprung gap and so evolving from the end 
of their main sequenc e to the base of t heir giant branch on 
a thermal timescale (jKolb et al.l ll997). As a consequence 
the mass-transfer rate should be too high for the standard 
soft X-ray tr a nsient explanation of the outbursts for which 
Ivan Paradiisl fl996) found that the average accretion rate 
on to the black hole in GRO J1655-40 should be less than 
1.26 x 10 -10 Mq yr -1 . Though we show that we can easily 
make models in which the secondary star is still on the main 
sequence, we still find that such a low mass-transfer rate is 
inco nsistent with any Roche- l obe fi lling system. 

iGreene. Bailvn fc Oroszl 1)20011 ) have carefully observed 
the secondary star in GRO J 1655-40 during quiescence. 
Their observations confirm that the secondary star is fill- 
ing its Roche lobe and are consistent with somewhat lower 
masses for both the secondary star and the black hole. In- 
deed we find that we must put the secondary star right at the 
lower end of their estimated mass range (Ma = 2.4±0.4 Mq) 
to fit its position in the Hertzspru ng-Russell diagram a t all 
well. The observed mass function (Shahbaz et al.|[l999h is 



f{M 2 ) = 



(Mi sini) 3 
(Mi + M 2 )'- 



= 2.73 ± 0.09 M , 



(35) 



wher e Mi is the mass of the blac k hole and i = 70.2° ± 
1.9° (|Greene. Bailvn fc Oroszll200 ll ) is the inclination of the 
binary orbit. In our models we choose the current secondary 
mass M2 and determine the current mass of the black hole 
Mi from the mass function and hence the total mass Mi + 
M2 of the system. We choose the initial secondary mass 
M21, the mass just before mass transfer begins. Assuming 
conservation of mass and angular momentum during mass 
transfer we have 



P(MiM 2 ) 3 = const 



(36) 



so we can calculate the initial period Pj needed to give the 
curre n t period of 2.62168 ± 0.00014 d (|van der Hooft et all 
1 19981 ) . IGreene. Bailvn fc Oroszl (|200ll ) measure a radius for 
the secondary star of 5.0 ± 0.3 Rq which we combine with 
their chosen effective temperature of 6,336 K to give a lumi- 
nosity L in the range 1.427 < log lo (Z//L ) < 1.667. 

We use the latest version of the Cambridge STARS code 
to construct detailed evolutionary models of the sec o ndary 
star. The code was originally written by lEggletonl l|l97ll. 
1972, 1973). The equation of state, which includes molecu- 
lar hydrogen, pressure ion i zation and coulomb interactions, 
is discussed bv lPols et al.l l|l995h . The initial composition is 
taken to be uniform with a hydrogen abundance X = 0.7, 
helium Y = 0.28 and metals Z = 0.02 with the m eteoritic 
mixture determined by lAnders fc Grevessd l|l989h . Hydro- 
gen burning is allowed by the p p chain and the CNO cycles . 
Reaction rates are taken from ICaughlan fc Fowleil (l988). 
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Figure 9. The evolutionary track in the Hertzsprung-Russell di- 
agram of a 2.5Mq star in a binary and which is forced to transfer 
mass to its companion when it fills its Roche lobe which occurs 
while crossing the Hertzsprung gap. The initial period and pri- 
mary mass are chosen so that the period and masses fit those 
of GRO J 1655-40 when the secondary reaches 2Mq. The dotted 
line and right-hand scale give the mass-transfer rate as a function 
of temperature. The point with error bars is the observed sec- 
ondary in GRO J1655-40 and the open circle is the model with 
M2 = 2 Mq at the observed period. 



Opaci ty ta bles are those calculated by llglesias fc Rogers! 
' l996t ) and I Alexander fc Ferguson] (| 19941 ). 



IFoellmi et al.l 1 2006h argue the source must be much 
nearer than typical estimates of 3.2 kpc. However, Bailyn 
(private communica tion) suggests that this is not possible. 
IFoellmi et alj (|2006l ) assumed that the secondary is a normal 
F6 IV star and so has the radius of an isolated F6 IV star. 
As we see below it is very easy to fit the secondary with a 
somewhat more massive and luminous but Roche-lobe filling 
star. 

After experimenting with a variety of initial conditions 
we find that, to fit its position in the Hertzsprung-Russell 
diagram, the current mass of the secondary star must be at 
the low end of the predicted range, between about 2.0 M 
and 2.1 Mq and that the mass of the black hole also lies 
at the low end of its predicted range, between 5.88 Mq and 
5.98 Mq. We illustrate two typical models in which the sec- 
ondary is in different evolutionary states. 

The first (Figsl9land ll0|) has a secondary star, of initial 
mass 2.5 Mq, and initial period 1.648 d. It fills its Roche lobe 
while evolving across the Hertzsprung gap. Evolution during 
this phase is on a thermal timescale so that the mass transfer 
rate is relatively high, 2.35 x 10 -7 Mq yr -1 , and lasts about 
3 x 10 6 yr with M > 10~ 7 M yr -1 . It fits the observations 
when the secondary's mass has fallen to 2 Mq and the black 
hole mass is 5.88 M . A total of AM = 0.5 M of material 
has passed through the disc. Similar results are found for 
lower masses but both the accretion rate and the time for 
which the black hole has been accreting fall. 



8 R. G. Martin, C. A. Tout & J. E. Pringle 





log 10 (T ef( /K) 

Figure 11. As Fig. [9] but with a secondary star initially of 
2.8 Mg. This time mass transfer begins on the main sequence 
when it is more gentle and lasts longer. 

The second (Figs [TT] and I12p has a secondary star, of 
initial mass 2.8 Mq, and initial period 1.481 d. This fills its 
Roche lobe while still evolving on the main sequence. Evo- 
lution at this phase is on a nuclear timescale so that the 
mass transfer rate is somewhat lower, 5.30 x 10 -9 Mq yr -1 , 
and lasts for about 5 x 10 7 yr. Once again this fits the ob- 
servations when M2 = 2Mq and M\ — 5.88 Mq. In this 
time AM = 0.8 Mq of material has passed through the 
disc. Secondary stars of initially higher mass transfer ma- 



terial at a higher rate but not for much longer. We also 
experimented with including convect ive overshooting (c.f. 
iRegos. Tout fc Wickramasinghelll998l ) to prolong the main- 
sequence lifetime but found very little difference in the mod- 
els that fit at these masses. 



7 DISCUSSION 

We have successfully modelled GRO J 1655-40 with mass 
transfer by Roche-lobe overflow. We need the masses of the 
black hole and its companion to be at the low end of their 
observed ranges. We find the average mass-transfer rate and 
the length of time that it lasts for different binary systems. 

We have described two typical, but not unique, mod- 
els in the previous section. If the companion is cross- 
ing the Hertzsprung gap a typical mass-transfer rate is 
2.35 x 10" 7 M Q yr" 1 . For this mass-transfer rate we predict 
an alignment timescale, equation (|34p . of i a ii gn = 7.8 x 10 yr. 
The accretion rate lasts for about 3 x 10 6 yr which is about 
3.8i a iign- For a secondary star still on the main sequence a 
typical mass-transfer rate is 5.30 x 10 -9 Mq yr -1 . The align- 
ment timescale for this lower rate is £ a ii g n = 2.5 x 10 7 yr. 
This accretion rate lasts for about 5 x 10 7 yr which is about 
2t a iign- We note that the alignment timescales also depend 
on the spin of the black hole which is not known. In these 
estimates we have taken a = 0.5. 

In calculating the timescales we have assumed conser- 
vative mass transfer and a constant accretion rate. Given 
that mass is expected to accumulate in the disc until an 
outburst the actual mass transfer might be expected to take 
place in bursts that last for less than the viscous timescale 
of the disc. During the intervening quiescence, when our 
equilibrium model can be applied, the mass-transfer rate 
would be lower and the alignment timescale correspondingly 
longer. However, a few tenths of a solar mass of material 
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must pass through the disc. The angular momentum of this 
material is much less than that of the black hole. Because 
faiign oc M~ 32/35 , taccrotc/taiign oc AM M~ 3/35 depends on 
M so changes in mass-transfer rate do not affect the ex- 
pected alignment very much. 

We can also entertain the possibility that the average 
ma ss-transfer rat e migh t be low enough to be consistent with 
the Ivan Paradiisl (1996) SXT model. In this case the align- 
ment timescale would be much longer. We do not have a 
consistent evolutionary model in this case but cannot expect 
the mass transfer to last longer than the main-sequence life- 
time of a 2 Mq star, about 8 x 10 8 yr. We do not however 
expect that such a low mass-transfer rate is possible from a 
Roche-lobe filling secondary star other than for a short time 
at the onset and the end of mass transfer. It is unlikely that 
we are catching the system in either of these brief periods 
and it is more likely that the condition for transience needs 
to be reconsidered. 

Alternatively the mass transfer might be highly non- 
conservative, though there is no observational evidence 
for this. In this case the alignment timescale would be 
much longer while the evolutionary timescales would be un- 
changed. Thus the models we have discussed would lead to 
the fastest alignment and this is just slow enough for us to 
expect to find the black hole still misaligned with the binary 
orbit. 

We see that, whether the secondary fills its Roche lobe 
on the main sequence or while evolving across the Hert- 
sprung gap, the accretion rate lasts for a few alignment 
timescales and so the fact that we see the black hole still 
misaligned with the orbital plane now is consistent. We are 
unable to distinguish whether the black hole is 15° from 
alignment or counteralignment. Given that the misalign- 
ment is likely to be created by a strong kick at the time 
of the supernova explosion we cannot expect either case to 
be more likely than the other at the start of the accretion. 

For constant mass transfer our model predicts exponen- 
tial decrease of the misalignment. We find that if r\ — 15° 
now and if mass transfer has been going on for 2t a ii gn then 
the misalignment would have initially been rjo = 148° . If the 
black hole has been accreting for 3.8 t a ii g n then the system 
would have initially been at 770 = 175°. However these are 
only indicative models and it is relatively easy to find others 
in which the accretion has been going on for less time, but 
still similar to the alignment timescale. We conclude that 
the system may or may not have been counter-aligned after 
its supernova kick but that it is most likely to be close to 
alignment rather than counteralignment now because oth- 
erwise the black hole would have had to form very close to 
7/0 — 180° and this is rather unlikely. 



8 CONCLUSIONS 

We have found that the secondary star is most likely evolv- 
ing either at the end of the main sequence or across the 
Hertzsprung gap. The mass-transfer rate is typically fifty 
times higher in the latter case but we find that, in both 
cases, the lifetime of the mass transfer state is at most a 
few times the alignment timescale. The fact that the black 
hole has not yet aligned with the orbital plane is therefore 
consistent with either model. The system may or may not 



have been counter-aligned after its supernova kick but we 
conclude that it is now more likely to be close to alignment 
rather than counteralignment. 
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